---
title: "OGGRAC"
output: 
  flexdashboard::flex_dashboard:
    storyboard: TRUE
    theme: 
      version: 4
      bg: "#101010"
      fg: "#FDF7F7" 
      primary: "#891812"
      navbar-bg: "#3ADAC6"
      base_font: 
        google: EB Garamond 
      heading_font:
        google: Sen
      code_font:
        google: 
          # arguments to sass::font_google() 
          family: JetBrains Mono
          local: false
    social: menu
    source_code: embed
---

```{r setup, include=FALSE}
library(flexdashboard)
library(sf)
library(leaflet)
library(plotly)
library(paletteer)
library(leaflet.extras)
library(leaflet.providers)
library(htmltools)
library(fontawesome)
library(plotly)
library(paletteer)
# Cambiar la fuente del título
tags$h1(style = "font-family: Arial, sans-serif;", "OGGRAC")

vul_muy_baja <- read_sf(dsn ="C:/Users/diego/Documents/OGRAC/MAPA_PRELIMINAR/wgs84/vulmuybaja.shp")
vul_baja <- read_sf(dsn ="C:/Users/diego/Documents/OGRAC/MAPA_PRELIMINAR/wgs84/vulbaja.shp")
vul_media <- read_sf(dsn ="C:/Users/diego/Documents/OGRAC/MAPA_PRELIMINAR/wgs84/vulmedia.shp")
vul_alta <- read_sf(dsn ="C:/Users/diego/Documents/OGRAC/MAPA_PRELIMINAR/wgs84/vulalta.shp")
vul_muy_alta <- read_sf(dsn ="C:/Users/diego/Documents/OGRAC/MAPA_PRELIMINAR/wgs84/vulmuyalta.shp")
sedes_unmsm <- readxl::read_excel("C:/Users/diego/Documents/OGRAC/sedes_unmsm/sedes_unmsm.xls")
dist_incendio_xy <- st_zm(read_sf(dsn ="C:/Users/diego/Documents/OGRAC/material_construccion/wgs84/riesgo_sismo.shp"), drop = TRUE)
sismos_sedes <- read_sf(dsn ="C:/Users/diego/Documents/OGRAC/MAPA_PRELIMINAR/PELIGROS_SISMOS_sedes.shp")
const_2d <- st_zm(read_sf(dsn ="C:/Users/diego/Documents/OGRAC/material_construccion/wgs84/mat_const_pol.shp"), drop = TRUE)
  
geomorfologia <- read_sf(dsn ="C:/Users/diego/Documents/OGRAC/fisico_espacial/wgs84/Geomorfologia_.shp")
geologia <- read_sf(dsn ="C:/Users/diego/Documents/OGRAC/fisico_espacial/wgs84/Geologia.shp")

labels1 <- sprintf(
  "<strong>%s</strong>",
  vul_muy_baja$EDIFICACIO
)%>% lapply(htmltools::HTML)
labels2 <- sprintf(
  "<strong>%s</strong>",
  vul_baja$EDIFICACIO
)%>% lapply(htmltools::HTML)
labels3 <- sprintf(
  "<strong>%s</strong>",
  vul_media$EDIFICACIO
)%>% lapply(htmltools::HTML)
labels4 <- sprintf(
  "<strong>%s</strong>",
  vul_alta$EDIFICACIO
)%>% lapply(htmltools::HTML)
labels5 <- sprintf(
  "<strong>%s</strong>",
  vul_muy_alta$EDIFICACIO
)%>% lapply(htmltools::HTML)





const_2d <- st_zm(read_sf(dsn ="C:/Users/diego/Documents/OGRAC/material_construccion/wgs84/mat_const_pol.shp"), drop = TRUE)

labels6 <- sprintf(
  "<b>Material Construcción:</b> %s<br><b>Número de pisos:</b> %s", 
  const_2d$MAT_PREDOM, 
  const_2d$PISOS_1
) %>% lapply(htmltools::HTML)

labels7 <- sprintf("<b>Número de pisos:</b>%s", const_2d$PISOS_1) %>% lapply(htmltools::HTML)
  pal1 <- colorFactor(
    palette = c("orange","red" ,"yellow"),
    levels = c( "Concreto armado", "Ladrillo", "Quincha")
  )
pal2 <- colorFactor(
    palette = c("#ffff01","#c5da03", "#619900" ,"#ffa901", "#ff2700" ),
    levels = c( "1", "2", "3", "4", "5")
  )

  




```



# TUTORIAL

## Column {data-width="300"}

### Información {data-height="300"}

En el siguiente podras acceder a los mapas realizados por el Equipo Técnico de la Oficina General de Gestión del Riesgo y Adaptación al Cambio Climático, asi como el "PLAN DE PREVENCIÓN Y REDUCCIÓN DEL RIESGO DE DESASTRES DE LA UNIVERSIDAD NACIONAL MAYOR DE SAN MARCOS 2024-2030".


<p style="text-align: center;">

<em><strong>[Oficina General de Gestión del Riesgo y Adaptación al Cambio Climático</strong>]{style="color: #FFFFFF;"}</em>

</p>

<p style="text-align: left;">

[<strong>GEOVISOR #1 <a href="https://diegomorantesuarez.github.io/OGGRAC.github.io/OGGRAC.html">https://diegomorantesuarez.github.io/OGGRAC.github.io/OGGRAC.html</a></strong>]{style="color: #FFFFFF;"}

</p>

<p style="text-align: left;">

[<strong>GEOVISOR #2 <a href="https://diegomorantesuarez.github.io/OGGRAC.github.io/OGGRAC2.html">https://diegomorantesuarez.github.io/OGGRAC.github.io/OGGRAC2.html</a></strong>]{style="color: #FFFFFF;"}

</p>

<p style="text-align: left;">

[<strong>DRIVE:
<a href="https://drive.google.com/drive/folders/1b3eeCcgKM2I6HGIGuGJFLKwZbYkvyhUe?usp=drive_link">https://drive.google.com/drive/folders/1b3eeCcgKM2I6HGIGuGJFLKwZbYkvyhUe?usp=drive_link</a></strong>]{style="color:  #FFFFFF;"}

</p>

<p style="text-align: left;">

 

</p>
## Column {data-width="500"}
### Partes del Geovisor {data-height="500"}

![](C:/Users/diego/Documents/OGRAC/introduccion/partes.png){.centered-image width="250"}

## Column {data-width="500"}

### Definición y Objetivos

![](C:/Users/diego/Documents/OGRAC/introduccion/definicion.png){.centered-image width="250"}

### Escala y temáticas

![](C:/Users/diego/Documents/OGRAC/introduccion/escala.png){.centered-image width="250"}



# SISMOS

## Column {data-width="500"}

### Vulnerabilidad Infraestructura

```{r}

m <- leaflet() %>%
  setView(lng = -77.082, lat = -12.056, zoom = 16)  %>%
  addTiles(
    urlTemplate = "https://server.arcgisonline.com/ArcGIS/rest/services/World_Imagery/MapServer/tile/{z}/{y}/{x}",
    attribution = 'Tiles &copy; Esri &mdash; Source: Esri, i-cubed, USDA, USGS, AEX, GeoEye, Getmapping, Aerogrid, IGN, IGP, UPR-EGP, and the GIS User Community',
    group = "Satelite"
  ) %>%
  addProviderTiles("OpenStreetMap", group = "Mapa" ) %>% 
  addPolygons(data = vul_muy_baja,  # Datos de los polígonos
              weight = 2,                   # Grosor del borde del polígono
              opacity = 1,                  # Opacidad del borde del polígono
              color = "green",        # Color del borde del polígono
              dashArray = "0.5",            # Tipo de línea del borde del polígono
              fillOpacity = 0.7,            # Opacidad del relleno del polígono
              popup = ~paste(
                "<h3 style='text-align:center;'>", vul_muy_baja$EDIFICACIO, "</h3><br/>",
                "<b>Valoración:</b>", vul_muy_baja$VALORACION, "<br/>",
                "<b>Material:</b>", vul_muy_baja$MATERIAL, "<br/>",
                "<b>Sistema Estructurado:</b>", vul_muy_baja$SIST_ESTRU, "<br/>",
                "<b>Intervención:</b>", vul_muy_baja$INTERVENCI, "<br/>",
                "<b>Estado de Conservación:</b>", vul_muy_baja$VALORA_CON, "<br/>",
                "<b>Año de Construcción:</b>", vul_muy_baja$AÑO_DE_CO
              ),  # Contenido del popup
              highlightOptions = highlightOptions(  # Opciones de resaltado
                weight = 2,                         # Grosor del borde al resaltar
                color = "#000",                     # Color del borde al resaltar
                dashArray = "",                     # Tipo de línea del borde al resaltar
                fillOpacity = 0.9,                  # Opacidad del relleno al resaltar
                bringToFront = TRUE),               # Traer al frente al resaltar
              label = labels1,                         # Etiquetas de los polígonos
              labelOptions = labelOptions(            # Opciones de las etiquetas
                style = list("font-weight" = "normal", padding = "3px 8px"),  # Estilo de la etiqueta
                textsize = "15px",                    # Tamaño del texto de la etiqueta
                direction = "auto"),                  # Dirección de la etiqueta
              group = "Descripción"                    # Agregar grupo para el control de capas
  ) %>% 
  addPolygons(data = vul_baja,  # Datos de los polígonos
              weight = 2,                   # Grosor del borde del polígono
              opacity = 1,                  # Opacidad del borde del polígono
              color = "chartreuse",        # Color del borde del polígono
              dashArray = "0.5",            # Tipo de línea del borde del polígono
              fillOpacity = 0.7,            # Opacidad del relleno del polígono
              popup = ~paste(
                "<h3 style='text-align:center;'>", vul_baja$EDIFICACIO, "</h3><br/>",
                "<b>Valoración:</b>", vul_baja$VALORACION, "<br/>",
                "<b>Material:</b>", vul_baja$MATERIAL, "<br/>",
                "<b>Sistema Estructurado:</b>", vul_baja$SIST_ESTRU, "<br/>",
                "<b>Intervención:</b>", vul_baja$INTERVENCI, "<br/>",
                "<b>Estado de Conservación:</b>", vul_baja$VALORA_CON, "<br/>",
                "<b>Año de Construcción:</b>", vul_baja$AÑO_DE_CO
              ),
              highlightOptions = highlightOptions(  # Opciones de resaltado
                weight = 2,                         # Grosor del borde al resaltar
                color = "#000",                     # Color del borde al resaltar
                dashArray = "",                     # Tipo de línea del borde al resaltar
                fillOpacity = 0.9,                  # Opacidad del relleno al resaltar
                bringToFront = TRUE),               # Traer al frente al resaltar
              label = labels2,                         # Etiquetas de los polígonos
              labelOptions = labelOptions(            # Opciones de las etiquetas
                style = list("font-weight" = "normal", padding = "3px 8px"),  # Estilo de la etiqueta
                textsize = "15px",                    # Tamaño del texto de la etiqueta
                direction = "auto"),                  # Dirección de la etiqueta
              group = "Descripción"                    # Agregar grupo para el control de capas
  ) %>% 
  addPolygons(data = vul_media,  # Datos de los polígonos
              weight = 2,                   # Grosor del borde del polígono
              opacity = 1,                  # Opacidad del borde del polígono
              color = "yellow",        # Color del borde del polígono
              dashArray = "0.5",            # Tipo de línea del borde del polígono
              fillOpacity = 0.7,            # Opacidad del relleno del polígono
              popup = ~paste(
                "<h3 style='text-align:center;'>", vul_media$EDIFICACIO, "</h3><br/>",
                "<b>Valoración:</b>", vul_media$VALORACION, "<br/>",
                "<b>Material:</b>", vul_media$MATERIAL, "<br/>",
                "<b>Sistema Estructurado:</b>", vul_media$SIST_ESTRU, "<br/>",
                "<b>Intervención:</b>", vul_media$INTERVENCI, "<br/>",
                "<b>Estado de Conservación:</b>", vul_media$VALORA_CON, "<br/>",
                "<b>Año de Construcción:</b>", vul_media$AÑO_DE_CO
              ),
              highlightOptions = highlightOptions(  # Opciones de resaltado
                weight = 2,                         # Grosor del borde al resaltar
                color = "#000",                     # Color del borde al resaltar
                dashArray = "",                     # Tipo de línea del borde al resaltar
                fillOpacity = 0.9,                  # Opacidad del relleno al resaltar
                bringToFront = TRUE),               # Traer al frente al resaltar
              label = labels3,                         # Etiquetas de los polígonos
              labelOptions = labelOptions(            # Opciones de las etiquetas
                style = list("font-weight" = "normal", padding = "3px 8px"),  # Estilo de la etiqueta
                textsize = "15px",                    # Tamaño del texto de la etiqueta
                direction = "auto"),                  # Dirección de la etiqueta
              group = "Descripción"                    # Agregar grupo para el control de capas
  ) %>% 
  addPolygons(data = vul_alta,  # Datos de los polígonos
              weight = 2,                   # Grosor del borde del polígono
              opacity = 1,                  # Opacidad del borde del polígono
              color = "orange",        # Color del borde del polígono
              dashArray = "0.5",            # Tipo de línea del borde del polígono
              fillOpacity = 0.7,            # Opacidad del relleno del polígono
              popup = ~paste(
                "<h3 style='text-align:center;'>", vul_alta$EDIFICACIO, "</h3><br/>",
                "<b>Valoración:</b>", vul_alta$VALORACION, "<br/>",
                "<b>Material:</b>", vul_alta$MATERIAL, "<br/>",
                "<b>Sistema Estructurado:</b>", vul_alta$SIST_ESTRU, "<br/>",
                "<b>Intervención:</b>", vul_alta$INTERVENCI, "<br/>",
                "<b>Estado de Conservación:</b>", vul_alta$VALORA_CON, "<br/>",
                "<b>Año de Construcción:</b>", vul_alta$AÑO_DE_CO
              ),
              highlightOptions = highlightOptions(  # Opciones de resaltado
                weight = 2,                         # Grosor del borde al resaltar
                color = "#000",                     # Color del borde al resaltar
                dashArray = "",                     # Tipo de línea del borde al resaltar
                fillOpacity = 0.9,                  # Opacidad del relleno al resaltar
                bringToFront = TRUE),               # Traer al frente al resaltar
              label = labels4,                         # Etiquetas de los polígonos
              labelOptions = labelOptions(            # Opciones de las etiquetas
                style = list("font-weight" = "normal", padding = "3px 8px"),  # Estilo de la etiqueta
                textsize = "15px",                    # Tamaño del texto de la etiqueta
                direction = "auto"),                  # Dirección de la etiqueta
              group = "Descripción"                    # Agregar grupo para el control de capas
  ) %>% 
  addPolygons(data = vul_muy_alta,  # Datos de los polígonos
              weight = 2,                   # Grosor del borde del polígono
              opacity = 1,                  # Opacidad del borde del polígono
              color = "red",        # Color del borde del polígono
              dashArray = "0.5",            # Tipo de línea del borde del polígono
              fillOpacity = 0.7,            # Opacidad del relleno del polígono
              popup = ~paste(
                "<h3 style='text-align:center;'>", vul_muy_alta$EDIFICACIO, "</h3><br/>",
                "<b>Valoración:</b>", vul_muy_alta$VALORACION, "<br/>",
                "<b>Material:</b>", vul_muy_alta$MATERIAL, "<br/>",
                "<b>Sistema Estructurado:</b>", vul_muy_alta$SIST_ESTRU, "<br/>",
                "<b>Intervención:</b>", vul_muy_alta$INTERVENCI, "<br/>",
                "<b>Estado de Conservación:</b>", vul_muy_alta$VALORA_CON, "<br/>",
                "<b>Año de Construcción:</b>", vul_muy_alta$AÑO_DE_CO
              ),
              highlightOptions = highlightOptions(  # Opciones de resaltado
                weight = 2,                         # Grosor del borde al resaltar
                color = "#000",                     # Color del borde al resaltar
                dashArray = "",                     # Tipo de línea del borde al resaltar
                fillOpacity = 0.9,                  # Opacidad del relleno al resaltar
                bringToFront = TRUE),               # Traer al frente al resaltar
              label = labels5,                         # Etiquetas de los polígonos
              labelOptions = labelOptions(            # Opciones de las etiquetas
                style = list("font-weight" = "normal", padding = "3px 8px"),  # Estilo de la etiqueta
                textsize = "15px",                    # Tamaño del texto de la etiqueta
                direction = "auto"),                  # Dirección de la etiqueta
              group = "Descripción"                    # Agregar grupo para el control de capas
  )%>%
  addPolygons(
    data = dist_incendio_xy[dist_incendio_xy$RIESGO == "ALTO", ],
    weight = 1,
    fillColor = "orange",
    color = "orange",
    fillOpacity = 1,
    group = "Riesgo"
  ) %>%
  addPolygons(
    data = dist_incendio_xy[dist_incendio_xy$RIESGO == "MEDIO", ],
    weight = 1,
    fillColor = "yellow",
    color = "yellow",
    fillOpacity = 1,
    group = "Riesgo"
  ) %>%
  addPolygons(
    data = dist_incendio_xy[dist_incendio_xy$NIV_VUL == "Vulnerabilidad Baja", ],
    weight = 1,
    fillColor = "chartreuse",
    color = "chartreuse",
    fillOpacity = 1,
    group = "Vulnerabilidad"
  ) %>%
  addPolygons(
    data = dist_incendio_xy[dist_incendio_xy$NIV_VUL == "Vulnerabilidad Media", ],
    weight = 1,
    fillColor = "yellow",
    color = "yellow",
    fillOpacity = 1,
    group = "Vulnerabilidad"
  )%>%
  addPolygons(
    data = dist_incendio_xy[dist_incendio_xy$NIV_VUL == "Vulnerabilidad Alta", ],
    weight = 1,
    fillColor = "orange",
    color = "orange",
    fillOpacity = 1,
    group = "Vulnerabilidad"
  )%>%
  addLayersControl(
    baseGroups = c("Descripción", "Vulnerabilidad", "Riesgo"),
    overlayGroups = c("Satelite", "Sedes"),  # Grupos de capas
    options = layersControlOptions(collapsed = F)
  ) %>%
  addLegend(
    position = "bottomright",
    colors = c("green","chartreuse", "yellow","orange","red"),
    labels = c("Muy baja","baja","Medio","alta","Muy alta"),
    title = "Vulnerabilidad Infraestructura"
  )%>%
  addLegend(
    position = "bottomright",
    colors = c( "yellow","orange"),
    labels = c("Medio","alta"),
    title = "Riesgo"
  ) %>%
  addControl(
    html = '<div style="position: absolute; bottom: 10px; left: 10px; background-color: rgba(32, 31, 96, 0.8); padding: 5px; border-radius: 5px;">
              <img src="http://oggrac.unmsm.edu.pe/wp-content/uploads/2021/02/Diseno-sin-titulo.png" 
                   width="200" height="64" style="background: transparent; border: none;">
            </div>',
    position = "bottomleft"
  ) %>%
    addMarkers(data = sedes_unmsm,
               lng = ~LONG,
               lat = ~LAT_1, 
               popup = ~as.character(Name), 
               label = ~as.character(Name),
               options = markerOptions(autoPan = FALSE),
               labelOptions = labelOptions(
                 noHide = F,  # Mantener la etiqueta visible al hacer clic fuera de ella
                 style = list("font-size" = "14px", "font-weight" = "bold")  # Estilo de la etiqueta
               ),
               icon = leaflet::makeIcon(
                 iconUrl = "C:/Users/diego/Documents/OGRAC/geovisorOGGRAC/sedes.png",  # Ruta a tu icono
                 iconWidth = 40, iconHeight = 40,  # Tamaño del icono
                 iconAnchorX = 15, iconAnchorY = 15,  # Anclaje del icono
                 popupAnchorX = 0, popupAnchorY = 0,  # Anclaje del popup
                 shadowUrl = NULL  # No se usa sombra en este caso
               ),
               group = "Sedes")
m
```

## Column {data-width="250"}

### Descripción

De acuerdo con lo resultante obtenido por el peligro y la vulnerabilidad, el nivel del escenario de riesgo de **"Nivel Alto"** está presente en seis sedes: San Fernando, Instituto Raúl Porras Barrenechea, Fundación - Museo Temple Radicati, Residencia Universitaria Julio Cesar Tello, Unidad de Postgrado de Derecho y Ciencia Política y Casona - Centro Cultural de San Marcos.

El riesgo de **"Nivel Medio"** se presenta en 7 sedes: Ciudad Universitaria, Sede San Juan de Lurigancho, Escuela Profesional de Tecnología Médica, Museo de Historia Natural y Escuela Profesional de Educación Física, Facultad de Medicina Veterinaria y Escuela Profesional de Nutrición y Escuela Profesional de Obstetricia.

Hay que considerar que las sedes de la **Casona - Centro Cultural San Marcos e Instituto Raúl Porras Barrenechea** son las que presentan un **mayor nivel de riesgo**. Las dos tienen más de 100 años de antigüedad. Casona - Centro Cultural San Marcos, actualmente se imparten talleres donde asisten niños y adolescentes (población vulnerable). Por lo que se recomienda que las clases sean reubicadas en otras sedes donde exista un mayor espacio, además se recomendaría hacer una pequeña capacitación de cómo trasladarse de manera ordenada en la presencia de un sismo.

La sede del **Instituto Raúl Porras Barrenechea** además de ser antigua, su segundo nivel tiene un estado de conservación malo donde actualmente se observa el cableado eléctrico expuesto y con presencia de humedad en la pared. Se recomienda realizar el saneamiento de paredes y renovar la geomembrana impermeable que protege al techo de la edificación, está presenta orificios donde se filtra el agua que daña las paredes.

Existe una **alta densidad de la comunidad universitaria** en la sede de Residencia Universitaria Julio Cesar Tello, San Fernando y la Escuela Profesional de Nutrición y Escuela Profesional de Obstetricia, lo que podría ocasionar un incremento en accidentes o pérdidas humanas al momento que ocurra un sismo. Por lo que se recomienda respetar o disminuir el aforo para evitar el hacinamiento en las instalaciones de las infraestructuras.

### Sismos en Sedes

```{r}
s <- leaflet() %>%
  setView(lng = -77.082, lat = -12.046, zoom = 11) %>%
  addTiles(
    urlTemplate = "https://server.arcgisonline.com/ArcGIS/rest/services/World_Imagery/MapServer/tile/{z}/{y}/{x}",
    attribution = 'Tiles &copy; Esri &mdash; Source: Esri, i-cubed, USDA, USGS, AEX, GeoEye, Getmapping, Aerogrid, IGN, IGP, UPR-EGP, and the GIS User Community'
  )  %>%
    addMarkers(data = sedes_unmsm,
               lng = ~LONG,
               lat = ~LAT_1, 
               popup = ~as.character(Name), 
               label = ~as.character(Name),
               options = markerOptions(autoPan = FALSE),
               labelOptions = labelOptions(
                 noHide = F,  # Mantener la etiqueta visible al hacer clic fuera de ella
                 style = list("font-size" = "14px", "font-weight" = "bold")  # Estilo de la etiqueta
               ),
               icon = leaflet::makeIcon(
                 iconUrl = "//upload.wikimedia.org/wikipedia/commons/thumb/3/3a/UNMSM_coatofarms_seal.svg/514px-UNMSM_coatofarms_seal.svg.png",  # Ruta a tu icono
                 iconWidth = 40, iconHeight = 40,  # Tamaño del icono
                 iconAnchorX = 15, iconAnchorY = 15,  # Anclaje del icono
                 popupAnchorX = 0, popupAnchorY = 0,  # Anclaje del popup
                 shadowUrl = NULL  # No se usa sombra en este caso
               ),
               group = "Sedes") %>% 
  addPolygons(data = sismos_sedes[sismos_sedes$Susc_Sismo == "Medio", ],
              weight = 1,
              fillColor = "yellow",
              color = "yellow",
              fillOpacity = 0.5,
              group = "Sismos") %>% 
  addPolygons(data = sismos_sedes[sismos_sedes$Susc_Sismo == "Alto", ],
              weight = 1,
              fillColor = "orange",
              color = "orange",
              fillOpacity = 0.5,
              group = "Sismos") %>% 
  addPolygons(data = sismos_sedes[sismos_sedes$Susc_Sismo == "Muy Alto", ],
              weight = 1,
              fillColor = "red",
              color = "red",
              fillOpacity = 0.5,
              group = "Sismos") %>% 
  addPolygons(data = sismos_sedes[sismos_sedes$Susc_Sismo == "FR", ],
              weight = 1,
              fillColor = "white",
              color = "white",
              fillOpacity = 0.5,
              group = "Sismos") %>% 
  addLegend(
    position = "bottomright",
    colors = c("white", "yellow","orange","red"),
    labels = c("Formación Rocosa","Medio","Alto","Muy Alto"),
    title = "Peligro Sismo", 
    opacity = 0.5
  )%>%
  addLayersControl(
    baseGroups = c("Sismos"),
    overlayGroups = c("Sedes"),  
    options = layersControlOptions(collapsed = T)
  )
s
```


# CARACTERÍSTICAS DE FÍSICO - ESPACIAL


## Columna 3 {data-width="400"}

### Estructura


```{r}
  

  
  # Crear el mapa Leaflet y agregar los polígonos
a <- leaflet() %>%
    setView(lng = -77.082, lat = -12.056, zoom = 16)  %>% addTiles(
    urlTemplate = "https://server.arcgisonline.com/ArcGIS/rest/services/World_Imagery/MapServer/tile/{z}/{y}/{x}",
    attribution = 'Tiles &copy; Esri &mdash; Source: Esri, i-cubed, USDA, USGS, AEX, GeoEye, Getmapping, Aerogrid, IGN, IGP, UPR-EGP, and the GIS User Community'
  ) %>%
    addPolygons(data = const_2d,  # Datos de los polígonos
                fillColor = pal1(const_2d$MAT_PREDOM),
                weight = 1,                   # Grosor del borde del polígono
                opacity = 1,                  # Opacidad del borde del polígono
                color = "black",        # Color del borde del polígono
                dashArray = "0.5",            # Tipo de línea del borde del polígono
                fillOpacity = 0.7,
                highlightOptions = highlightOptions(  # Opciones de resaltado
                  weight = 2,                         # Grosor del borde al resaltar
                  color = "#000",                     # Color del borde al resaltar
                  dashArray = "",                     # Tipo de línea del borde al resaltar
                  fillOpacity = 0.9,                  # Opacidad del relleno al resaltar
                  bringToFront = TRUE),               # Traer al frente al resaltar
                label = labels6,                         # Etiquetas de los polígonos
                labelOptions = labelOptions(            # Opciones de las etiquetas
                  style = list("font-weight" = "normal", padding = "3px 8px"),  # Estilo de la etiqueta
                  textsize = "15px",                    # Tamaño del texto de la etiqueta
                  direction = "auto"),                  # Dirección de la etiqueta
                group = "Mat. Construcción"                    # Agregar grupo para el control de capas
    )%>%
    addPolygons(data = const_2d,  # Datos de los polígonos
                fillColor = pal2(const_2d$PISOS_1),
                weight = 1,                   # Grosor del borde del polígono
                opacity = 1,                  # Opacidad del borde del polígono
                color = "black",        # Color del borde del polígono
                dashArray = "0.5",            # Tipo de línea del borde del polígono
                fillOpacity = 0.7,
                highlightOptions = highlightOptions(  # Opciones de resaltado
                  weight = 2,                         # Grosor del borde al resaltar
                  color = "#000",                     # Color del borde al resaltar
                  dashArray = "",                     # Tipo de línea del borde al resaltar
                  fillOpacity = 0.9,                  # Opacidad del relleno al resaltar
                  bringToFront = TRUE),               # Traer al frente al resaltar
                label = labels6,                         # Etiquetas de los polígonos
                labelOptions = labelOptions(            # Opciones de las etiquetas
                  style = list("font-weight" = "normal", padding = "3px 8px"),  # Estilo de la etiqueta
                  textsize = "15px",                    # Tamaño del texto de la etiqueta
                  direction = "auto"),                  # Dirección de la etiqueta
                group = "Cantidad Pisos"                    # Agregar grupo para el control de capas
    ) %>%
    addMarkers(data = sedes_unmsm,
               lng = ~LONG,
               lat = ~LAT_1, 
               popup = ~as.character(Name), 
               label = ~as.character(Name),
               options = markerOptions(autoPan = FALSE),
               labelOptions = labelOptions(
                 noHide = F,  # Mantener la etiqueta visible al hacer clic fuera de ella
                 style = list("font-size" = "14px", "font-weight" = "bold")  # Estilo de la etiqueta
               ),
               icon = leaflet::makeIcon(
                 iconUrl = "//upload.wikimedia.org/wikipedia/commons/thumb/3/3a/UNMSM_coatofarms_seal.svg/514px-UNMSM_coatofarms_seal.svg.png",  # Ruta a tu icono
                 iconWidth = 40, iconHeight = 40,  # Tamaño del icono
                 iconAnchorX = 15, iconAnchorY = 15,  # Anclaje del icono
                 popupAnchorX = 0, popupAnchorY = 0,  # Anclaje del popup
                 shadowUrl = NULL  # No se usa sombra en este caso
               ),
               group = "Sedes") %>%
    addLayersControl(
      baseGroups = c("Cantidad Pisos", "Mat. Construcción"),
      overlayGroups = c("Sedes"),  # Grupos de capas
      options = layersControlOptions(collapsed = F)
    )%>% 
  addLegend(
    position = "bottomright",
    colors = c("#ffff01","#c5da03", "#619900" ,"#ffa901", "#ff2700" ),
    labels = c( "1 piso", "2 pisos", "3 pisos", "4 pisos", "5 pisos"),
    title = "Cantidad de Pisos", 
    opacity = 0.5) %>% 
  addLegend(
    position = "bottomright",
    colors = c("orange","red" ,"yellow"),
    labels = c("Concreto armado", "Ladrillo", "Quincha"),
    title = "Material Construcción", 
    opacity = 0.5)%>%
  addControl(
    html = '<div style="position: absolute; bottom: 10px; left: 10px; background-color: rgba(32, 31, 96, 0.8); padding: 5px; border-radius: 5px;">
              <img src="http://oggrac.unmsm.edu.pe/wp-content/uploads/2021/02/Diseno-sin-titulo.png" 
                   width="200" height="64" style="background: transparent; border: none;">
            </div>',
    position = "bottomleft"
  )
a
```







## Columna 3 {data-width="400"}

### Geomorfología y Geología

```{r}

geomorfologia_label <- sprintf(
  "<strong>%s</strong>",
  geomorfologia$etiqueta
)%>% lapply(htmltools::HTML)


  # Crear el mapa Leaflet y agregar los polígonos
a <- leaflet() %>%
    setView(lng = -77.020, lat = -12.056, zoom = 12)  %>% addTiles(
    urlTemplate = "https://server.arcgisonline.com/ArcGIS/rest/services/World_Imagery/MapServer/tile/{z}/{y}/{x}",
    attribution = 'Tiles &copy; Esri &mdash; Source: Esri, i-cubed, USDA, USGS, AEX, GeoEye, Getmapping, Aerogrid, IGN, IGP, UPR-EGP, and the GIS User Community',
    group = "Satélite"
  ) %>%
  addProviderTiles("OpenStreetMap", group = "Mapa" ) %>%
    addMarkers(data = sedes_unmsm,
               lng = ~LONG,
               lat = ~LAT_1, 
               popup = ~as.character(Name), 
               label = ~as.character(Name),
               options = markerOptions(autoPan = FALSE),
               labelOptions = labelOptions(
                 noHide = F,  # Mantener la etiqueta visible al hacer clic fuera de ella
                 style = list("font-size" = "14px", "font-weight" = "bold")  # Estilo de la etiqueta
               ),
               icon = leaflet::makeIcon(
                 iconUrl = "C:/Users/diego/Documents/OGRAC/geovisorOGGRAC/sedes.png",  # Ruta a tu icono
                 iconWidth = 40, iconHeight = 40,  # Tamaño del icono
                 iconAnchorX = 15, iconAnchorY = 15,  # Anclaje del icono
                 popupAnchorX = 0, popupAnchorY = 0,  # Anclaje del popup
                 shadowUrl = NULL  # No se usa sombra en este caso
               ),
               group = "Sedes")%>% 
  addPolygons(data = geomorfologia[geomorfologia$etiqueta == "Río", ],
              weight = 1,
              fillColor = "#73dfff",
              color = "#73dfff",
              fillOpacity = 0.5,
              group = "Geomorfología") %>% 
  addPolygons(data = geomorfologia[geomorfologia$etiqueta == "P-at", ],
              weight = 1,
              fillColor = "#eb9beb",
              color = "#eb9beb",
              fillOpacity = 0.5,
              group = "Geomorfología")%>% 
  addPolygons(data = geomorfologia[geomorfologia$etiqueta == "Pl-al", ],
              weight = 1,
              fillColor = "#9eaad7",
              color = "#9eaad7",
              fillOpacity = 0.5,
              group = "Geomorfología")%>% 
  addPolygons(data = geomorfologia[geomorfologia$etiqueta == "RCL-ri", ],
              weight = 1,
              fillColor = "#456e8a",
              color = "#456e8a",
              fillOpacity = 0.5,
              group = "Geomorfología")%>% 
  addPolygons(data = geomorfologia[geomorfologia$etiqueta == "RCL-rs", ],
              weight = 1,
              fillColor = "#d5de52",
              color = "#d5de52",
              fillOpacity = 0.5,
              group = "Geomorfología")%>% 
  addPolygons(data = geomorfologia[geomorfologia$etiqueta == "RCL-rvs", ],
              weight = 1,
              fillColor = "#a6803f",
              color = "#a6803f",
              fillOpacity = 0.5,
              group = "Geomorfología")%>% 
  addPolygons(data = geomorfologia[geomorfologia$etiqueta == "RM-ri", ],
              weight = 1,
              fillColor = "#4d4dd6",
              color = "#4d4dd6",
              fillOpacity = 0.5,  
              group = "Geomorfología")%>% 
  addPolygons(data = geomorfologia[geomorfologia$etiqueta == "RM-rvs", ],
              weight = 1,
              fillColor = "#a83800",
              color = "#a83800",
              fillOpacity = 0.5,
              group = "Geomorfología")%>% 
  addPolygons(data = geomorfologia[geomorfologia$etiqueta == "T-m", ],
              weight = 1,
              fillColor = "#90d487",
              color = "#90d487",
              fillOpacity = 0.5,
              group = "Geomorfología")%>% 
  addPolygons(data = geologia[geologia$UNIDAD == "Depósito aluvial", ],
              weight = 1,
              fillColor = "#c29ed7",
              color = "#c29ed7",
              fillOpacity = 0.5,
              group = "Geología")%>% 
  addPolygons(data = geologia[geologia$UNIDAD == "Depósito aluvial, fluvial", ],
              weight = 1,
              fillColor = "#3ba149",
              color = "#3ba149",
              fillOpacity = 0.5,
              group = "Geología")%>% 
  addPolygons(data = geologia[geologia$UNIDAD == "Depósito coluvial", ],
              weight = 1,
              fillColor = "#ebdc71",
              color = "#ebdc71",
              fillOpacity = 0.5,
              group = "Geología")%>% 
  addPolygons(data = geologia[geologia$UNIDAD == "Depósito fluvial", ],
              weight = 1,
              fillColor = "#73b2ff",
              color = "#73b2ff",
              fillOpacity = 0.5,
              group = "Geología")%>% 
  addPolygons(data = geologia[geologia$UNIDAD == "Depósito marino", ],
              weight = 1,
              fillColor = "#857036",
              color = "#857036",
              fillOpacity = 0.5,
              group = "Geología")%>% 
  addPolygons(data = geologia[geologia$UNIDAD == "Formación Atocongo", ],
              weight = 1,
              fillColor = "#ffebaf",
              color = "#ffebaf",
              fillOpacity = 0.5,
              group = "Geología")%>% 
  addPolygons(data = geologia[geologia$UNIDAD == "Formación Pamplona", ],
              weight = 1,
              fillColor = "#996559",
              color = "#996559",
              fillOpacity = 0.5,
              group = "Geología")%>% 
  addPolygons(data = geologia[geologia$UNIDAD == "Grupo Casma - Formación Chilca", ],
              weight = 1,
              fillColor = "#7971eb",
              color = "#7971eb",
              fillOpacity = 0.5,
              group = "Geología")%>% 
  addPolygons(data = geologia[geologia$UNIDAD == "Grupo Casma - Formación Quilmaná", ],
              weight = 1,
              fillColor = "#65f783",
              color = "#65f783",
              fillOpacity = 0.5,
              group = "Geología")%>% 
  addPolygons(data = geologia[geologia$UNIDAD == "Grupo Morro Solar - Formación Herradura", ],
              weight = 1,
              fillColor = "#8c679c",
              color = "#8c679c",
              fillOpacity = 0.5,
              group = "Geología")%>% 
  addPolygons(data = geologia[geologia$UNIDAD == "Grupo Morro Solar - Formación Marcavilca", ],
              weight = 1,
              fillColor = "#f79681",
              color = "#f79681",
              fillOpacity = 0.5,
              group = "Geología")%>% 
  addPolygons(data = geologia[geologia$UNIDAD == "Grupo Puente Piedra - Formación Ancón", ],
              weight = 1,
              fillColor = "#4ecca0",
              color = "#4ecca0",
              fillOpacity = 0.5,
              group = "Geología")%>% 
  addPolygons(data = geologia[geologia$UNIDAD == "Grupo Puente Piedra - Formación Puente Inga", ],
              weight = 1,
              fillColor = "#c4f2a2",
              color = "#c4f2a2",
              fillOpacity = 0.5,
              group = "Geología")%>% 
  addPolygons(data = geologia[geologia$UNIDAD == "Grupo Puente Piedra - Formación Ventanilla", ],
              weight = 1,
              fillColor = "#365196",
              color = "#365196",
              fillOpacity = 0.5,
              group = "Geología")%>% 
  addPolygons(data = geologia[geologia$UNIDAD == "Sin denominación", ],
              weight = 1,
              fillColor = "#f79cdf",
              color = "#f79cdf",
              fillOpacity = 0.5,
              group = "Geología")%>% 
  addPolygons(data = geologia[geologia$UNIDAD == "Súper Unidad Patap", ],
              weight = 1,
              fillColor = "#a3405c",
              color = "#a3405c",
              fillOpacity = 0.5,
              group = "Geología")%>% 
  addPolygons(data = geologia[geologia$UNIDAD == "Súper Unidad Santa Rosa", ],
              weight = 1,
              fillColor = "#4d8074",
              color = "#4d8074",
              fillOpacity = 0.5,
              group = "Geología") %>%
    addLayersControl(
      baseGroups = c( "Geomorfología","Geología"),
      overlayGroups = c("Sedes", "Satélite"),  # Grupos de capas
      options = layersControlOptions(collapsed = F)
    )
a
```


## Columna 3 {.tabset .tabset-fade data-width="300"}

### Descripción

**Número de pisos de la edificación**

Según el Método Meseri las edificaciones con mayor altura o mayor número de pisos tienen mayor probabilidad de propagación y es más difícil su extinción.

**Material de construcción de la edificación**

El Método de Meseri indica que las estructuras más resistentes al fuego, reducen el nivel de propagación del fuego en las edificaciones Según este método, las construcciones de hormigón y similares tienen alta resistencia, los metales protegidos mediana resistencia y las construcciones de madera tienen una baja resistencia. Así observamos que cerca del 75% de las edificaciones son de concreto armado, mientras que de ladrillo son 11% y las de madera el 7%, mientras que las áreas no construidas corresponden al 7% de las edificaciones, con un riesgo nulo.

**Área de la construcción de la edificación**

Las edificaciones mientras más extensas sean, requieren de una mayor exigencia de resistencia al fuego; según el DBSI 2009 el riesgo es directamente proporcional al área construida o al volumen construido. Por lo que el riesgo se incrementa mientras más extensa sea el área construida de la edificación. De acuerdo con con *el tercer gráfico* la UNMSM en su ámbito de estudio tiene 7 áreas entre sin construir y construidas hasta 100 metros cuadrados, tiene también 21 edificaciones con un área menor a 500 metros cuadrados, 20 edificaciones que se encuentran entre 500 y 1000 metros cuadrados, 48 edificaciones que están entre 1000 y 5000 metros cuadrados y 3 con un área mayor a los 5000 metros cuadrados.


### Grafico

```{r}


df <- data.frame(
  Número.de.pisos = c("Sin construcción", "1-2 pisos", "3 pisos", "4 pisos", "5 pisos a más"),
  Nivel.de.Propagación = c("Nulo, Sin construcción", "Bajo, 1-2 pisos", "Medio, 3 pisos", "Alto, 4 pisos", "Muy alto, 5 pisos a más"),
  Cantidad.de.edificaciones = c(7 , 43 , 26 , 20 , 3 ),
  Porcentaje = c(7.07, 43.44, 26.26, 20.2, 3.03)
)



df2 <- data.frame(
  Material.de.Construcción = c("Área sin construcción", "Concreto armado", "Ladrillo", "Adobe", "Madera, Quincha y otros"),
  Nivel.de.Resistencia = c("Nula, Área sin construcción", "Muy alta, Concreto armado", "Alta, Ladrillo", "Media, Adobe", "Baja, Madera, Quincha y otros"),
  Cantidad.de.edificaciones = c(7 , 74 , 11 , 0 , 7 ),
  Porcentaje2 = c(7.07, 74.75, 11.11, 0, 7.07)
)

df3 <- data.frame(
  Área.de.Construcción = c("0-100", "101-500", "501-1000", "1001-5000", "Mayor a 5000"),
  Nivel.de.exigencia = c("Muy bajo, 0-100 ", "Bajo, 101-500", "Medio, 501-1000", "Alto, 1001-5000", "Muy alto, Mayor a 5000"),
  Cantidad.de.edificaciones = c(7 , 21 , 20 , 48 , 3 ),
  Porcentaje3 = c(7.07, 21.21, 20.2, 48.49, 3.03)
)



# Asignar colores
colores <- paletteer_c("ggthemes::Orange-Gold", 5) 
colores2 <- paletteer_c("ggthemes::Green-Gold", 5)
colores3 <- paletteer_c("ggthemes::Blue-Teal", 5)


fig <- plot_ly() %>% 
  add_pie(data = df, labels = ~Nivel.de.Propagación, values = ~Cantidad.de.edificaciones ,name = "N° Pisos" , type = 'pie', marker = list(colors = colores), domain = list(row = 0, column = 0)) %>% 
  add_pie(data = df2, labels = ~Nivel.de.Resistencia, values = ~Cantidad.de.edificaciones ,name = "Mat. Construcción", type = 'pie', marker = list(colors = colores2), domain = list(row = 1, column = 0))%>% 
  layout(title = "Características constructivas", showlegend = F, 
         grid=list(rows=2, columns=1),
         xaxis = list(showgrid = FALSE, zeroline = FALSE, showticklabels = FALSE),
         yaxis = list(showgrid = FALSE, zeroline = FALSE, showticklabels = FALSE)) %>% layout(
           paper_bgcolor = '#101010',  # Color de fondo del papel (fondo exterior)
           plot_bgcolor = 'white',    # Color de fondo del gráfico (fondo interior)
           title = list(
             text = "Características constructivas",
             font = list(color = "white")  # Color del título
           ),
           font = list(color = "white")
         )
fig

```

### **Geomorfología y Geología** 

|             **Geomorfología**              | **Descripción**                             |
|:-----------------------------------:|:----------------------------------|
| ![](images/1.png){width="30" height="19"}  | Depósito aluvial                            |
| ![](images/2.png){width="30" height="19"}  | Depósito aluvial, fluvial                   |
| ![](images/3.png){width="30" height="19"}  | Depósito coluvial                           |
| ![](images/4.png){width="30" height="19"}  | Depósito fluvial                            |
| ![](images/5.png){width="30" height="19"}  | Depósito marino                             |
| ![](images/6.png){width="30" height="19"}  | Formación Atocongo                          |
| ![](images/7.png){width="30" height="19"}  | Formación Pamplona                          |
| ![](images/8.png){width="30" height="19"}  | Grupo Casma - Formación Chilca              |
| ![](images/9.png){width="30" height="19"}  | Grupo Casma - Formación Quilmaná            |
| ![](images/10.png){width="30" height="19"} | Grupo Morro Solar - Formación Herradura     |
| ![](images/11.png){width="30" height="19"} | Grupo Morro Solar - Formación Marcavilca    |
| ![](images/12.png){width="30" height="19"} | Grupo Puente Piedra - Formación Ancón       |
| ![](images/13.png){width="30" height="19"} | Grupo Puente Piedra - Formación Puente Inga |
| ![](images/14.png){width="30" height="19"} | Grupo Puente Piedra - Formación Ventanilla  |
| ![](images/15.png){width="30" height="19"} | Sin denominación                            |
| ![](images/16.png){width="30" height="19"} | Súper Unidad Patap                          |
| ![](images/17.png){width="30" height="19"} | Súper Unidad Santa Rosa                     |

|                **Geología**                | **Descripción**                                          |
|:-----------------------------------:|:----------------------------------|
| ![](images/g1.png){width="30" height="19"} | Río                                                      |
| ![](images/g2.png){width="30" height="19"} | P-at (Vertiente o piedemonte aluvio - torrencial)        |
| ![](images/g3.png){width="30" height="19"} | Pl-al (Llanura o planicie aluvial)                       |
| ![](images/g4.png){width="30" height="19"} | RCL-ri (Colina y lomada en roca intrusiva)               |
| ![](images/g5.png){width="30" height="19"} | RCL-rs (Colina y lomada en roca sedimentaria)            |
| ![](images/g6.png){width="30" height="19"} | RCL-rvs (Colina y lomada en roca volcano - sedimentaria) |
| ![](images/g7.png){width="30" height="19"} | RM-ri (Montaña en roca intrusiva)                        |
| ![](images/g8.png){width="30" height="19"} | RM-rvs (Montaña en roca volcano sedimentaria)            |
| ![](images/g9.png){width="30" height="19"} |  T-m (Terraza marina)                                    |

```{r}

```